Dynamic networks and directed percolation (28 January 2009) 

Roni Parshani, 1 Mark Dickison, 2 Reuven Cohen, 3 H. Eugene Stanley, 2 and Shlomo Havlin 1 

1 Minerva Center & Department of Physics, 
Bar-Ilan University, Ramat Gan, Israel 
2 Center for Polymer Studies, Boston University, 
Boston, Massachusetts 02215, USA 
3 Department of Mathematics, Bar-Ilan University, Ramat Gan, Israel 
(Dated: 28 January 2009 — pdcsh.tex) 

Abstract 

We introduce a model for dynamic networks, where the links or the strengths of the links change 
over time. We solve the model by mapping dynamic networks to the problem of directed percola- 
tion, where the direction corresponds to the evolution of the network in time. We show that the 
dynamic network undergoes a percolation phase transition at a critical concentration p c , which 
decreases with the rate r at which the network links are changed. The behavior near criticality is 
universal and independent of r. We find fundamental network laws are changed, (i) For Erdos- 
Renyi networks we find that the size of the giant component at criticality scales with the network 
size N for all values of r, rather than as iV 2 / 3 . (ii) In the presence of a broad distribution of dis- 
order, the optimal path length between two nodes in a dynamic network scales as iV 1 / 2 , compared 
to iV 1 / 3 in a static network. 
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Network theory ha 8 answered many q nestio» 8 eoncerning S tatie networks Q fl Q, fl 
[5( but many real networks are dynamic in the sense that their links, or the strengths of 
their links, change with time. For example, in social networks friendships are formed and 
dissolved, while in communication networks, such as the Internet, the load (weight) on the 
links changes continually. The challenges posed by such dynamic networks are beginning 
to be addressed. For example, Kempe et al. have studied algorithms for broadcasting or 
gossiping in dynamic networks (| while Volz and Meyers have studied the epidemic SIR 
model on dynamic networks 3]. 

Fundamental questions that have been extensively studied in static networks are still 
open for dynamic networks. Here we ask: (i) Is there a critical concentration of links for 
which the dynamic network undergoes a percolation phase transition, above which order 
N of the network nodes are still connected and below which the network breaks into small 
clusters? (ii) If so, what is the percolation threshold for which the transition occurs, and 
how does it depend on the dynamics? (iii) What are the properties near criticality? 

Consider an iV-node network with M links that change over time. Each link has a lifetime 
r drawn from a Poisson distribution, and is replaced by a new link between two randomly- 
selected nodes after its lifetime r expires (Fig. [1]). We define a unit time step as the time 
required for a walker to traverse one link; we also assume that when traversing the network 
one can only remain for a limited time at each node. Even if there is no path between nodes 
A and B at a specific time, a walker traversing the network may be able to pass from point 
A to point B because new links are continually appearing. Likewise, even if a path between 
A and B exists at a given time it may be disconnected before a walker is able to traverse it 

(Fig. up. 

In order to represent the time evolution of the network we add to our network represen- 
tation another axis, which corresponds to a time axis. Each layer along the time dimension 
corresponds to a network state at some time t [Fig. [2(a)]. Our new network representation 
includes all the information regarding the network dynamics. 

We now argue that percolation on a dynamic network is equivalent to the problem of 



directed percolation (DP) in infinite dimensions 10|, llll]. A percolation process can be 
described as an invasion of a liquid in a porous media where the flow between two neighboring 
pores is permitted with some probability]? and blocked with probability (1— p). This process 
undergoes a phase transition at the critical probability p = p c . For p > p c the liquid will 
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FIG. 1: A five-node dynamic network at different time steps. The weights W are defined in the 
text. 

spread throughout the media, while for p < p c the liquid will be constrained to a finite 
area. At exactly p c the liquid fills the incipient infinite cluster (called the giant component 
in network theory). Directed percolation (DP), a special kind of percolation, is a process in 



which the spread is limited to a direction that is defined to be the longitudinal axis 



While non-directed percolation has been successfully 
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161 ]. the analog of DP in networks 



is widely used for studying network stability 
has not been previously studied. 

In a dynamic network the direction has the following interpretation. The evolution of the 
network in time plays the role of the direction and the additional time axis presented in our 
model acts as an additional longitudinal (vertical) axis that is the hallmark of DP [Fig. EJ^a) 
and Fig. [T]. The correspondence between percolation on dynamic networks and DP on 
regular networks allows us to apply the results known from DP at criticality to dynamic 
networks. 

Networks can be regarded as infinite dimensional structures since no spatial constraints 
exist. Therefore we expect the critical properties of dynamic networks to be the same as 
DP in infinite dimensions. The relevant critical properties for DP are 10, |jjj]: (i) D(t), the 
number of nodes reached at time t, scales as D(t) ~ t; S(t) (ii) The giant component size, 
scales as S(t) ~ t 2 ; and (iii) P s (t), the survivability (the probability of reaching layer t when 
growing a cluster), scales as P s (t) ~ £ _1 . Figures [2](b) and[3]^a) present simulation results 
confirming these scaling relations. 

To learn about the size- dependent properties of dynamic networks we determine the DP 
properties as a function of the network size N, rather than as a function of t. In DP at 
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FIG. 2: (Color online) (a) The y axis corresponds to the time dimension and each successive 
layer of two rows corresponds to a network configuration at a different time step (c.f. Fig.l.) 
The full red line path is the optimal path (sum of its weights is minimal) between node A and 
node E, even though a shorter path exists, shown as a black dashed line, (b) Simulation results 
for the survivability P s (t) and its cutoff, at criticality, for dynamic networks of different network 
sizes. From left to right: N = 100, 400, 1600, 6400, 25600. The data collapse shown in the inset 
demonstrates that P s (t) in dynamic networks is universal when scaled by N 1 ^ 2 . 

criticality, the infinite dimensional relationship between w, the width in the transverse axes, 
and t, the length in the longitudinal axes, is w ~ t 1 ! 2 . The upper critical dimension d c 
is the lowest dimension for which the system has the properties of an infinite dimensional 
system. For DP this value is d c = 4 + 1 (1 corresponds to the longitudinal axis), so the 
relation between the system size at the upper critical dimension and the size of a dynamic 
network is given by N ~ w 4 (the power 4 comes from the 4 transverse dimensions of d c ). 
Since w ~ t 1 ! 2 we conclude that: 

t ~ N 1/2 . (1) 

Therefore for a dynamic network of size N at criticality, P s {t) decays exponentially after a 
time t x , with t x ~ N 1 ' 2 Q- Fi gure El^b) presents simulation results for the survivability 
of the giant cluster in a dynamic network at criticality. The figure shows that, for different 
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values of N and t > t x , P s (t) ~ as expected from DP in infinite dimensions. The 
exponential decay for t x > N 1 ^ 2 can also be seen, in agreement with Eq. ([T|). The inset 
of Fig. Wfo) shows the collapse of survivability data after scaling by iV 1 / 2 , supporting again 
Eq. p. 

The size of the giant component as a function of N is derived by substituting Eq. ([T]) in 
the DP relation S(t) ~ t 2 . 

S(N) ~ N. (2) 

Figure [3]^b) presents simulation results illustrating this scaling relationship, as well as a 
corresponding relationship for static networks, where S(N) is known to scale as S(N) ~ N 2 / 3 
|s 5 2]- Ps(t) for static networks is also known to decay exponentially after a time t x ~ iV 1 / 3 . 
The two systems clearly have different behavior and properties at criticality, and thus belong 
to two different universality classes. 

Next we show that the behavior of dynamic networks at criticality is universal and inde- 
pendent of the rate in which the links are changed [inset of Fig. H](a)]. We do find, however, 
that the critical concentration, p c , for which the phase transition occurs depends on the av- 
erage lifetime of the links [Fig. IH(a)]. The formula for p c on Erdos-Renyi dynamic networks 
is given by 

p c =(k)-^((k) + c(k)) 1 -^ (3) 

where (r) is the average rate with which the links change and (k) is the average degree. 
Fig. H](b) presents simulation results for S(N) at different values of p indicating the network 
undergoes a phase transition at some critical value of p. The inset of Fig. IH(a) shows that 
for several different values of r S(N) ~ iV at p c , as expected due to the universality of the 
critical exponent in dynamic networks. The agreement of our simulations with the formula 
for p c [Eq. (j3J)] is shown in Fig. H](a). 

The correspondence to DP can also predict the general scaling of the optimal path in 
a dynamic network with a broad distribution of disorder. In a network where weights are 
assigned to links, the optimal path between any two nodes is defined as the path along which 
the sum of the weights is minimal. In the limit of a broad distribution of disorder, Ref. [3] 
has shown that, at criticality, the optimal path exists mainly along the giant cluster. There- 
fore for static networks the optimal path length scales with the average distance between 
nodes on the percolation cluster: £ opt ~ N 1 ^ 3 . In our dynamic network model the average 
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FIG. 3: (a) The cluster size S(t) scales as t 2 . (b) Simulation results (dots) supporting the relation 
S(N) ~ N (Eq. 3) (upper dashed line), the cluster size in dynamic networks, compared to the 
known S(N) ~ iV 2 / 3 (bottom dashed line), the cluster size in static networks. 

distance between nodes on the percolation cluster scales as (£) ~ iV 1 / 2 , suggesting that in 
dynamic ER networks the optimal path scales as 

L P t ~ N 1 ' 2 . (4) 

Figure [5] shows simulation results for the optimal path length in a dynamic network compared 
to a static network. The results for dynamic networks are in full agreement with Eq. (jlj). 

What makes the results in a dynamic network so different from the static case? The 
difference lies in the number of available configurations. While in static networks the per- 
colation cluster is composed from paths built from N network nodes, in dynamic networks 
the network is represented by N 3 / 2 nodes {^(J. The evolution of the network over time 
generates many more possible configurations, enabling the percolation cluster to become 
much larger. Substituting N' = N 3 / 2 in the percolation cluster formula for static networks, 
S{N) ~ iV 2 / 3 , yields S(N') ~ S{N 3 / 2 ) ~ N, which further confirms our results for S{N) in 
dynamic networks. 

The explanation for the optimal path being longer in directed networks is that it can 
include loops. For example, an optimal path reaching some node A at t — t\ may find it 




FIG. 4: (a) Simulation fit to the formula for p c as a function of the rate (r) at which the links 
change. The value of p c was calculated by finding the value of p for which S(N) gives a straight 
line. In the inset we show the value of S(N) at p c for several different values of (r). In all cases the 
slope equals 1 as predicted by the critical exponent for dynamic networks, (b) Simulation results 
for S(N) in a network with {k) = 3 and (r) = 2/3 are presented for different values of p showing 
that the network undergoes a phase transition at p c . 

necessary to return to that node at t2 > ti if at £2 node A is more optimally connected to 
the destination node. In a static network the optimal path does not include loops since any 
link connecting node A at £2 was also available at t\ therefore any loop will only increase 
the total weight. 

Representing a dynamic network as a directed network [Fig. [2^ a)] composed of N 3 ^ 2 nodes 
allows the "same" node to be counted more then once in the percolation cluster, therefore 
requiring that the distinct number of nodes on the percolation cluster also scale with N. To 
determine the number of different nodes of the original network in a component of size M 
on the directed network, consider the following argument: The links between consecutive 
layers of the directed network are chosen randomly. Therefore, each link leads to a random 
node in the original network independently and with uniform distribution. The probability 
to reach a new node by following a link, assuming that D nodes have already been visited, 
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FIG. 5: The optimal path for strong disorder scales as £(N) ~ N 1 / 2 in dynamic networks compared 
to £(N) ~ N 1 / 3 in static networks. 

is 1 — D/N. The expected number of distinct nodes E(D) reached after A links have been 
followed from the starting node is therefore E(D X )/N = E(D X -i) + E(l — D X -i/N). This 
reduces to E(D X ) = 1 + (1 - l/N)E(D X -i) which indicates that for large M 

ffi-l = l-(l-^) 1 '» !-«-/«. (5) 

Thus, when the size of a component in the directed network is of order iV a finite fraction 
of the visited nodes are new and the size of the induced component on the original network 
is also of order N. 

In summary, we introduced a model for dynamic networks which was solved by a compar- 
ison with directed percolation in infinite dimensions. The DP longitudinal axis is mapped to 
the time axis along which the dynamic network evolves. We showed that dynamic networks 
exhibit different properties and critical exponents near criticality. Therefore they belong to 
a different universality class than static networks. While in static networks S(N), the size of 
the giant component at criticality, scales as S(N) ~ iV 2//3 , in dynamic networks S(N) ~ N. 
Even though the properties of dynamic networks are universal and independent of the rate r 
at which the links are changed, the critical concentration, p c , for which the phase transition 
occurs depends on r. We also showed that the optimal path in dynamic networks scales as 
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^opt ~ iV 1 / 2 , compared to £ opt ~ A^ 1 / 3 in static networks. 

We thank the Israel Science Foundation, Yeshaya Horowitz Association and The Center 
for Complexity Science, EU Daphnet project, DOE and ONR for financial support, and N. 
Madar and R. Bartsch for discussions. 



[1] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 (2002). 

[2] R. Albert et al., Nature 406, 378 (2000). 

[3] R. Cohen et al., Phys. Rev. Lett. 85, 4626 (2000). 

[4] R. Pastor-Satorras and A. Vespignani Phys. Rev. Lett. 86, 3200 (2001). 
[5] M. E. J. Newman, Phys. Rev. E 66, 016128 (2002). 

[6] D. Kempe, J. Kleinberg, and A. Kumar, Proc. 32nd ACM STOC. (2000); Proc. 33rd ACM 
STOC. (2001). 

[7] E. Volz and L. A. Meyers, Proc. Royal Soc. B (2007); J. Royal Soc. Interface, 

doi:10.1098/rsif.2008.0218 (2008). 
[8] P. Erdos and A. Renyi, Publ. Math. (Debrecen), 6, 290 (1959); A. Publ. Math. Inst. Hung. 

Acad. Sci. 5, 1760 (1960). 
[9] B. Bollobas, Random Graphs (Academic Press, Orlando, 1985). 
[10] W. Kinzel, in Percolation Structures and Processes, edited by G. Deutscher, R. Zallen, and J. 

Adler (Hilger, Bristol, 1983). 
[11] H. Hinrichsen, Adv. Phys. 49, 815 (2000). 
[12] D. S. Callaway et al., Phys. Rev. Lett. 85, 5468 (2000). 
[13] R. Cohen et al., Phys. Rev. E 66, 036113 (2002). 
[14] N. Schwartz et al., Phys. Rev. E 66, 015104 (2002). 
[15] E. Lopez et al., Phys. Rev. Lett. 99, 188701 (2007). 
[16] R. Cohen et al., Phys. Rev. Lett. 86, 3682 (2001). 

[17] Reversing the relation from DP we get t(w) ~ w 2 , implying that if we limit the distance w 
(instead of the conventional way in DP of limiting the time t) we will obtain a cutoff in the 
known critical exponents after time t x ~ w 2 . The equivalent for networks would be a cutoff 
after a crossover time t x ~ N 1 / 2 (based on Eq. ([1])). 

9 



[18] For Erdos-Renyi static networks the percolation threshold is given by p c = (k)^ 1 . Since at 
criticality the average number of new nodes reachable from a given node approaches 1, p(k) = 
1. However, for percolation on our dynamic network model where a node can be revisited 
since the links may change, there are two limits. First, for networks having links with a long 
average lifetime, a node reached in a percolation process is connecetd to an additional link 
(through which it was reached) on top of the (k) links to which an average node is connected. 
Such a node has a higher effective average degree that depends on the probability c((k)), to 
return to a previous node. The second limit is for networks changing constantly. In this case 
the degree of a node reached in a percolation process is a random variable with no memory, 
and its effective average degree remains (k). Therefore, the value of p c in dynamic networks 
depends on the rate at which the links change and includes a combination of both limits. 
A good approximation for the probability of a node visited in a percolation process to be 
revisited is obtained when considering the first order contribution of more than one path to 
return to the original node after two steps. This leads to c((k)) ~ (1 — p 2 )( k \ The theoretical 
line in Fig. 0] was calculated according to this approximation. 

[19] L. A. Braunstein et al., Phys. Rev. Lett. 91, 168701 (2003). 

[20] The directed network representation of a dynamic network as shown in Fig. [21(a) includes tN 
nodes where t x ~ N 1 / 2 according to Eq. ([2]). 



10 



